# Replication Package for Labor Market Effects of Nuclear Power Plants
# Authors: Duha T. Altindag (dza0013@auburn.edu), Reem El Cheikh Taha (rme0017@auburn.edu), Jennifer U. Jones (juj2102@columbia.edu), and R. Alan Seals, Jr. (ras0029@auburn.edu)
# This replication package contains all the necessary files and instructions to replicate the results presented in the paper and appendix.

## Contents
- `do_files/`: Directory containing the Stata script that replicates the results in the paper
   - 'ReplicationMain_Figures.do': Stata script that replicates all the figures in the paper
   - 'ReplicationMain_Tables.do': Stata script that replicates all the tables in the paper
   - 'ReplicationAppendix.do': Stata script that replicates all the figures and tables in Appendices A & B
   - 'SC weights builder.do': Stata script that produces the weights in Appendix Table B3 (Note that completion of this code may take several hours).

- `clean_data/`: Directory containing the datasets used in the study.
	-'clean_data/Final Panel.dta': Panel dataset of all U.S.A. counties including economic outcomes and nuclear power plants details from data sources described in the text, covering the period 1969-2019
	-'clean_data/Considered Panel.dta': Panel dataset of the counties that built their first reactor between 1974 & 1978, including economic outcomes and nuclear power plants details from data sources described in the text, covering the period 1969-2019
	-'clean_data/Considered Panel.dta': same as above, but on commuting zone level instead of county level as described in the text
	-'clean_data/Considered Panel_gov finance.dta': Panel dataset of the counties that built their first reactor between 1974 & 1978, including government finance variables from data sources described in the text, covering the period 1969-2019
	-'clean_data/Final SC panel (Emp).dta': Panel dataset of the treated counties and their synthetic controls, where the outcome to be matched is employment-to-population
	-'clean_data/Final SC panel (Wages).dta': Panel dataset of the treated counties and their synthetic controls, where the outcome to be matched is the log of per capita wages and salaries 
	-'clean_data/power stations.dta': hand-collected dataset on all nuclear reactors in the U.S.A., covering major events such as the year of announcement, construction, operation, and more as described in the text
	-'clean_data/proxies.dta': dataset on proxies for the location characteristics published by the NRC as described in the text.
	-'clean_data/panel of neighbors of both treated and control according to population.dta': Panel dataset of the neighboring counties (according to the highest population) of both treated and control, as described in the text
	-'clean_data/cb_2019_us_county_20m.dta': U.S.A. state and county coordinates to plot the map in Figure 3.
	-'clean_data/Leave one out': Panel dataset of the treated counties and their synthetic controls, leaving out one control county at a time to produce Appendix Figures B5, B6, B7, & B8
	-'clean_data/falsification test_Donor treated': Panel dataset of donor counties and their synthetic controls, assuming they are treated to produce Appendix Figures B5, B6, B7, & B8
	-'clean_data/07819-0004-Data-card_image.txt': cbp poll as described in the text
	-'clean_data/Synthetic Control: folders including the weights of the SC of each treated county where employment to population or log of wages per capita are the outcomes

- `results/`: Directory to save the output.

- `temp_data/`: Directory to save temporary data.

- `README.txt`: This file provides an overview and instructions for replication.

## Requirements
- Software: Stata 16.1
- Required user-written Stata Packages: 
	-esttab: Author: Ben Jann (2007). estout: Stata module to make regression tables. Statistical Software Components S456501, Boston College Department of Economics. https://ideas.repec.org/c/boc/bocode/s456501.html
	-psmatch2: Leuven, E., & Sianesi, B. (2003). PSMATCH2: Stata module to perform full Mahalanobis and propensity score matching, common support graphing, and covariate imbalance testing. Statistical Software Components S432001, Boston College Department of Economics.
	-reghdfe: Correia, Sergio. 2014. "REGHDFE: Stata module to perform linear or instrumental-variable regression absorbing any number of high-dimensional fixed effects." Statistical Software Components S457874, Boston College Department of Economics. 
	-ftools: Correia, S. (2016). FTOOLS: Stata module to provide alternatives to common Stata commands optimized for large datasets. Statistical Software Components S458213, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s458213.html
 	-spmap: Pisati, M. (2007). SPMAP: Stata module to visualize spatial data. Statistical Software Components S456812, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s456812.html
 	-eventstudyinteract: Sun, Liyang. 2021. EVENTSTUDYINTERACT: Stata module to implement the interaction weighted estimator for an event study. Statistical Software Components S458978, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s458978.html 
 	-avar: Baum, C. F., & Schaffer, M. E. (2013). AVAR: Stata module to perform asymptotic covariance estimation for iid and non-iid data robust to heteroskedasticity, autocorrelation, 1- and 2-way clustering, and common cross-panel autocorrelated disturbances. Statistical Software Components S457689, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s457689.html 
 	-coefplot: Jann, B. (2013). COEFPLOT: Stata module to plot regression coefficients and other results. Statistical Software Components S457686, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s457686.html
 	-did_multiplegt_dyn:de Chaisemartin, C., Ciccia, D., D'Haultfoeuille, X., Knau, F., MalÃ©zieux, M., & Sow, D. (2023). DID_MULTIPLEGT_DYN: Stata module to estimate event-study Difference-in-Difference (DID) estimators in designs with multiple groups and periods, with a potentially non-binary treatment that may increase or decrease multiple times. Statistical Software Components S459222, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s459222.html
 	-regsave: Reif, J. (2008). REGSAVE: Stata module to save regression results to a Stata-formatted dataset. Statistical Software Components S456964, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s456964.html
 	-synth (version 0.0.7): Abadie, A., Diamond, A., & Hainmueller, J. (2011). SYNTH: Stata module to implement Synthetic Control Methods for Comparative Case Studies. Statistical Software Components S457334, Boston College Department of Economics. Available at: https://ideas.repec.org/c/boc/bocode/s457334.html
 	-gtools: Caceres Bravo, M. (2018). GTOOLS: Stata module to provide a fast implementation of common group commands. Statistical Software Components S458514, Boston College Department of Economics. 


## Instructions

1.	Unzip the compressed files in the “stata files/clean_data” 
2.	Open the 'do_files/ReplicationMain_Figures.do' & 'do_files/ReplicationMain_Tables.do' & 'do_files/ReplicationAppendix.do' & 'do_files/SC weights builder.do' files in Stata.
3.	Change the working directory to the directory in which you are storing this data and code package.
4. 	Install the above-mentioned Stata packages. 
5.	Run 'do_files/ReplicationMain_Figures.do' 'do_files/ReplicationMain_Tables.do' and 'do_files/ReplicationAppendix.do' to run the analysis script in the main text and Appendix.
6.      Note that the code ('do_files/SC weights builder.do') that produces the weights in Appendix Table B3 is commented out in the 'do_files/ReplicationAppendix.do' (line 1179). It can take several hours to complete.



